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2 Abstract 

^ We consider the UqSU{2) invariant spin-| XX Z quantum spin chain at roots of unity 

O q = exp(^^^), corresponding to different minimal models of conformal field theory. 

Q We conduct a numerical investigation of correlation functions of UqSU{2) scalar two- 

point operators in order to find which operators in the minimal models they correspond 
to. Using graphical representations of the Temperley-Lieb algebra we are able to deal 
with chains of up to 28 sites. Depending on q the correlation functions show different 



I characteristics and finite size behaviour. For m = |, which corresponds to the Lee- Yang 

edge singularity, we find the surface and bulk critical exponent — |. Together with the 
O known result in the case m = 3 (Ising model) this indicates that in the continuum limit 

the two-point operators involve conformal fields of spin ^i^^- For other roots of unity q 
the chains are too short to determine surface and bulk critical exponents. 
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1. Introduction 

We consider two-point correlation functions for a class of one-dimensional quantum 
models on a chain of N sites defined in terms of the hamiltonian [f , 2] 

JV-l 

H = -J2e,. (f.f) 

8 = 1 

Here the , i = f , 2, . . . , — f , are generators of a Temperley-Lieb algebra rjv(<?) 
defined by the relations [3] 

e, = {q + q~^) e, {1.2a) 
e,±i e, = e, {1.2b) 
e, ej = ej e, {j ^ i ± 1) . {1.2c) 

In this article we numerically compute correlation functions of two-point operators that 
are elements of this algebra. Representing the generators by 

e. = -^(<^Hi + ^M+i + ^^(^M+i-1) + ^^(^f-^Hi; 

these two-point operators become UqSU{2) invariant generalizations of the scalar 
operators Here erf, erf, and erf are Pauli matrices acting on site i. In the 

representation (1.3) H turns into a UqSU{2) symmetric spin-| XX Z Heisenberg 
hamiltonian. Throughout this article N is even to ensure that the ground state is 
unique. 

If q is real, the model is in a massive phase and in the limit N ^ oo correlation 
functions of local operators can be computed [4]. For the SU{2) symmetric model 
{q = 1) the correlation function of the operator crlcrf has been computed numerically 
using exact diagonalization [5, 6, 7] or Bethe ansatz [8]. The same correlation function 
has been calculated analytically in [9] for the periodic XX Z chain (where the UqSU{2) 
symmetry is lost). For this case numerical results can be found in [10]. 

In this work we consider the hamiltonian H and the two-point scalar operators 
gk,i (as defined in section 2) in terms of Temperley-Lieb generators. This includes the 
UqSU{2) symmetric case (1.3). Since in our calculations we use the Temperley-Lieb 
relations only, the results apply also to other representations besides this spin-| chain 
[11]. We take q = exp(^^^) a root of unity. In this case the model is massless and the 
spectrum of H contains the spectrum of hamiltonians corresponding to conformal field 
theories with central charge [12, 13, 14] 

c = l- . • (1.4) 

Although in the representation (1.3) the hamiltonian is not hermitian, its spectrum is 
real [1]. 
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The two-point scalar operators are non-local and their continuum limit is not known 
for general m. However, for m = 2 and m = 3 the correlation functions have been derived 
analytically using quotients of the Temperley-Lieb algebra [2, 15]. In the former case 
the correlation functions vanish. In the latter case H can be related to the hamiltonian 
of the Ising model and one finds two different non- vanishing correlation functions, whose 
continuum limit is given by correlation functions of operators with conformal dimensions 
{h,h) = (|,0) and (0,|). For m = 5 there is a quotient of TN{q) which gives a three 
state Potts model, but the correlation function has not been computed in this case. 
In this article we present a numerical investigation of the correlation functions in the 
cases m = 3, 5, 1 corresponding to Ising model, three state Potts model, and Lee- Yang 
edge singularity respectively. The first case is included to estimate the accuracy of 
the numerical results. For m = 5 we find four different correlation functions {gk,i)^ 
depending on k and / odd or even. In the case m = | there is only one correlation 
function. We attempt to find critical exponents and to identify conformal fields that 
correspond to the continuum limit of the two-point scalar operators. Due to a symmetry 
of the two-point scalar operator and from the known results in the Ising case one expects 
conformal fields with spin This can be confirmed in the case m = |. 

For the computation of the correlation functions we use graphical representations 
of rjv(<?) on a path space and on boundary diagrams [16, 11]. In the representation 
on boundary diagrams one can easily restrict the configuration space to a space, which 
is related to the UqSU{2) scalars. For integer m we can use a corresponding path 
representation space with a dimension that is further reduced. Whereas in the spin 
space approach (1.3) a reduction of the representation to the space of UqSU{2) scalars 
cannot be achieved easily. For the representation on boundary diagrams we can handle 
chains of up to 24 sites (independent of m), and up to 28 sites with the path space 
representation for m = 5. We extrapolate the data for — )• oo to correlations in the 
semi-infinite geometry, that is we keep one end of the chains fixed. From that we try to 
compute bulk and surface exponents. 

The article organizes as follows. In section 2 we define the two-point scalar operator 
and review some results from [2]. The graphical representations of the Temperley-Lieb 
algebra we use are introduced in section 3. In section 4 we describe the calculation of 
ground state and correlation functions. We analyze the results of the computation for 
the cases m = 3, 5, | in section 5, numerical data for the semi-infinite chain are listed 
in an appendix. Section 6 contains our conclusions. 
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2. Two-point scalar operators 

Two-point UqSU{2) scalar operators have been derived in [2]. They can be defined 
recursively in terms of Temperley-Lieb generators e^: 

gk,k+i = ek - (q + q-')-' l<k<N-l (2.1a) 

gk,i = -qgk,ngn,i - q~^ gn,i gk,n l<k<n<l<N {2.1b) 

gi,k =q-^gk,i l<k<l<N. (2.1c) 

Note that these operators are non-local. Interchanging q and one can define another 
set of such operators, but they have the same expectation values. For q = I and the 
representation (1.3) the two-point scalar operator reduces to the SU{2) scalar operator 

glT^ = -\^kOi . (2.2) 

For q = exp(^^^) a root of unity quotients of the Temperley-Lieb algebra can be 
used to relate the general hamiltonian (1.1) to other quantum spin chains. Specifically, 
for q = exp(^) the Temperley-Lieb algebra can be represented such that H becomes 
the hamiltonian of the Ising model with = y sites. With this representation the two- 
point scalar operators can be expressed as products of fermionic operators and their 
ground state expectation values can be computed explicitly [2]: 

{g2j,2k) = (2.3a) 
(52,-i,2fc-i) = (2.36) 

{gU2k-r) =-^^E sm(vr^^,)cos(vr^^(^--)) (2.3c) 

= ^^Esm(vr^^^)cos(vr^^(,--)).(2.3^) 

These expressions reduce in the limit Z — )• 00 to 

\/2 1 1 

{gk,i) = {-J T - J—t) k even, / odd (2.4a) 

1 1 

{gk,i) = ^ ij^ + j^) k odd, / even . (2.46) 

We will use these results to estimate errors of the numerical computations. 

For q = exp(^) there exists a different quotient of the Temperley-Lieb algebra such 
that the hamiltonian H turns into the one of a self-dual three state Potts quantum chain 
with Y sites and free boundary conditions [2]. Using this representation the operators 
gk,i can be rewritten in terms of two different local parafermionic operators, which 
correspond to respectively even or odd sites of the original quantum chain (1.1). The 
hamiltonian and expressions for the g^^i are explicitly given in [2]. The parafermionic 
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operators are already known from [17]. However, an analytic solution of the correlation 
functions using the parafermions has not been achieved. 

3. Two representations of the Temperley— Lieb algebra 

For the numerical computations we use two different graphical representations of the 
Temperley-Lieb algebra. We describe them in the following and explain how they can 
be employed to calculate correlation functions. 

3.1. Path representation 

First let q be generic. Then one can define an action of the Temperley-Lieb algebra on 
a vector space S with an orthonormal basis of vectors Vk that are labelled by + 1 
numbers k = (A;o, . . . , A;jv), subject to the conditions 

k > 

ko = kN = (3.1) 
k = k_i±l/2 i = l,...,N. 

On this vector space a generator of rjv(<?) can be represented by [11, 16] 



with 

k' = {ki,...,k,_i,kl,k,+i,...,kN) ■ (3.3) 
Here we use the definition of the q-numbei 



Defining 



q-q ^ 



r'''= " ]- (34) 



S has the dimension Fq^^ [1]. 

The vectors can be interpreted as paths of a Bratteli diagram [1, 16]. This diagram 
describes the fusion process of the N spin-| representations of UqSU{2) attached to 
each site of the spin chain. For the Bratteli diagram ki gives the spin of an irreducible 
representation of UqSU{2) that appears when one decomposes the tensor product of 
the spin-A;i_i multiplet attached to the first i — 1 sites with the doublet of site i. Thus 
k^ gives the total spin, and the condition k^ = selects the UqSU{2) scalars. In this 
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context the UqSU{2) symmetric generators (1.3) are found to act on the path space 
according to (3.2). 

Next we consider the case q = exp(^^^) = exp(^) with r and s coprime integers. 
Since [s]q = 0, in this case a basis of the path representation space of rjv(<?) is given by 
the vectors Vk with 

We denote this vector space as Sq C S. This reduction is reflected by the UqSU{2) 
representations for q a root of unity [1]. In this case indecomposable but reducible 
representations appear in the decomposition of the spin configuration space. Only the 
spin zero representations whose paths in the Bratteli diagram are restricted according 
to (3.5) remain as irreducible representations. Their number, i.e. the dimension of iSg, 
is given by [1] 

0(N) _ -p(iV) -p(iV) , -p(iV) -p(iV) , -p(iV) -p(iV) , /o fiX 

— -"-0 J- (m + l)-l "T J- m + 1 J- 2(m + l)-l + J- 2(m + l) ^ 3{m + l)-l ~^ " ' V^-^) 

and thus depends on q = exp(^^^). Since the ground state of the hamiltonian (1.1) is 
known to be non-degenerate, we can restrict the representation space for the numerical 
calculation of the correlation functions to S respectively Sq. The action of hamiltonian 

(1.1) and two-point scalar operator (2.1a-c) on the path space follow from the definition 

(3.2) . 

3.2. Representation on boundary diagrams 

In [2] a different diagrammatic approach to the calculation of correlation functions {gk,i) 
has been proposed, that uses the regular representation of the Temperley-Lieb algebra 
on boundary diagrams [16]. Boundary diagrams are given by N non-crossing lines that 
connect N upper and N lower points such that each point is connected to one other 
point. For the generators we have the definition: 



I 8 + 1 



(3.7) 



The diagrams corresponding to words of the algebra are obtained as follows. For the 
composition of two words the corresponding diagrams are stacked on top of each other 
and the lines ending in the lower points of the first diagram are joined with the lines 
ending in the corresponding upper points of the second diagram. Any closed line 
appearing in this process is discarded from the diagram and replaced by a factor q^q~^ . 
This refiects the Temperley-Lieb relation (1.2a). The other Temperley-Lieb relations 
(1.26) and (1.2c) can easily be verified by drawing the corresponding diagrams. In 
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this way elements of the Temperley-Lieb algebra can be represented faithfully as linear 
combinations of boundary diagrams [16]. 

Let us restrict this regular representation of rjv(<?) to the one-sided ideal X generated 
by the word 6163 .. . ejv-i- We have 

KJ KJ KJ 

6163 . . . ejv_i = • • • . (3.8) 

The ideal X is known to represent the space of UqSU{2) scalars for generic q [16]. Hence 
the ground state can be represented by a vector in X. The dimension of X (as of S) is 
given by Fq^^ The ideal is realized by all diagrams with disconnected upper and lower 
parts where the lower part is the one in (3.8) and the upper part has any configuration 
of y non-intersecting lines between the N upper points. 

To calculate expectation values one needs a scalar product on X. Define a transposed 
diagram as the upside-down refiected diagram. With this definition the scalar product 
(si|s2) of two arbitrary elements Si and S2 of X is defined [2]: 

WW W 

Sl^S2 = {si\s2) ■ ■ ■ . (3.9) 

That is, the scalar product is given by the factor that one obtains removing any loop 
from the diagrams which represent Si"^S2- 

With the definition of boundary diagrams given above and the scalar product 
(3.9), the expectation values {v\gk^i\v) of the correlation operators can be computed 
by multiplying diagrams, that encode the Temperley-Lieb relations. This way it can 
easily be seen that {v\gk^i\v) is real for any vector v in X . 

For q a root of unity quotients of the Temperley-Lieb algebra have been given 
in [15]. The calculation of the ground state average of g^^i through evaluation of the 
corresponding diagrams according to (3.9) is independent of the quotienting, since it 
depends on the Temperley-Lieb relations (1.2a-c) only. In principle the quotients can be 
used to eliminate certain vectors from X. However, this is not useful for the numerical 
calculations, since the elimination procedure would have to be carried out after each 
application of a generator to a vector of X. 

4. Computational method 

Let us describe the numerical calculation of the ground state expectation values {gk,i) of 
the two-point operators (2.1a-c). We employ the representations of section 3 so that we 
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do not have to refer to the spin configuration space. In this article we present the cases 
q = exp(^^^) with m = 3, m = 5, and m = |, but the computation can be performed 
analogously for any m. 

First we generate a basis of the configuration space X (boundary diagrams) or Sq 
(path representation) and the matrix of the hamilton operator in the respective basis. 
In general the path representation has the advantage of a smaller representation space 
and an orthonormal basis of this space. Note that the matrices are not hermitian 
(although the eigenvalues of H are real) so that we cannot use the Lanczos algorithm. 
We employ a power method to find the ground state: repeated application of the matrix 
to an appropriate start vector projects onto the eigenstate with the eigenvalue of highest 
absolute value. It requires about 2000 iterations and gives the eigenenergy with a relative 
error of less than 10~^^ (10~^® for small N). We have only applied this method in the 
cases where H is represented by a real matrix. The complex case requires more storage 
capacity and multiplications. For the representation on boundary diagrams the matrix 
realizing H is real independent of the value of q. With the path representation we have 
a real hamiltonian for integer m. For the non-unitary model given by m = | the path 
representation gives a complex hamiltonian so that we use only the boundary diagrams. 
Note that for m = | the norm of the ground state is negative according to (3.9). 

For m = 3 and m = 5 the ground state energy was checked against Bethe ansatz 
calculations [13, 18]. Our results (obtained up to A^=28) are identical to the results 
listed in [13]. Also, using the formulas given in [12, 13], extrapolation from the finite size 
energies gives the correct values (1.4) of the conformal charge. We find c = 0.50(1) for 
m = 3, c = 0.79(1) for m = 5, and c = —4.4(1) for m = | using the BST extrapolation 
algorithm [19, 20]. 

We find that the ground state has positive parity, i.e. the following property: The 
contribution of a specific vector Vk from Sq (a word from X) to the ground state is the 
same as the one of the vector labelled by k' = (A;jv, ^Af-i, • • • , ^i) (the word with left and 
right ends exchanged). 

For the computation of the correlation functions we make explicit use of the recursion 
in (2.1a-c). We also use that the expectation values {v\gk^i\v) of the operators g^^i have 
the symmetry 

{v\9k,i\^) = {^\9N+i-i,N+i-k\v) (4.1) 

for any vector v with positive or negative parity. This can most easily be seen employing 
the boundary diagrams as a regular representation and using the reality of {v\gk^i\v). 

On a workstation we are able to calculate the correlations for chains of up to = 24 
sites with the representation using boundary diagrams (independent of q). In this 
case the representation space has dimension Fq^^^ = 208012. The limit is set by the 
time required for the computation of the scalar products (3.9). For the path space 
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representation the storage capacity limits the calculations. In this case the dimension 
of the representation space depends on m. We can handle = 40 sites for m = 3 
(dim 5, = = 524288) and iV = 28 sites for m = 5 {n^^% = 797162). 

From the known results in the Ising case (m = 3) we can determine the numerical 
errors of the computed correlations. The results are exact to 14 digits for = 14 and to 
11 digits for N = 28. The two approaches give the same values up to 12 digits [N = 24). 

We use the BST algorithm to extrapolate the correlation functions for finite N to 
the values of the semi-infinite chain [N — )• oo), i.e. we extrapolate the values {gk,i) for 
different N with k and / fixed. With the given precision of our numerical calculation we 
need at least the results of 6 different chain lengths to obtain reliable results (4 digits) 
with the BST extrapolation. This can be seen from the convergence behaviour of the 
BST algorithm. In the case m = 3 we can also check the extrapolated results using 
(2.4a-6) and find the same accuracy. Hence for the semi-infinite chain we can give {gk,i) 
with A;, / < 18 (m = 5) and A;, / < 14 (m = 2/3). These values are listed in the appendix. 



5. Results 



Here we discuss the correlation functions, i.e. the ground state expectation values {gk,i)-, 
for the three cases q = exp(^^^) with m = 3, 5, |. We consider the extrapolated results 
for the semi-infinite chain. For q on the unit circle the model is massless and thus we 
expect long-range order with power-law decay of the correlation functions. In the bulk 
limit (without boundary infiuence) we expect 

(^m) oc (5.1) 

with the standard definition of a bulk critical exponent x. From this and the relation 
gk,i = gi,k for k < I we find that 
m — 1 



m -\- L 

The continuum limit of the correlation functions of the semi-infinite chain 
corresponds to two-point functions of a conformal field theory in the half plane with 
central charge given by (1.4). In these minimal models the highest weights are found to 
be [21] 

[(m -|- l)p — mqY — 1 
4:m[m -\- 1) 

for q and p integers. Comparing (5.2) with (5.3) we are forced to take 

X = /ii,3(m) . (5-4) 

This implies that in the continuum limit the correlation function {gk,i) is realized by 
two-point functions of linear combinations of conformal fields with dimension (/ii,3, 0) 



hp,q{m) = , , , , (5.3) 
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or (0,/ii^3). For m = 3 the analytic result is known and one actually finds Majorana 
fields with spin | [2, 22]. Nonetheless we describe the analysis of the numerical results 
for this simple case, in order to test methods of finding information about the continuum 
limit of the correlation functions. By means of the numerical investigations we try to 
confirm (5.4) for other values of m. 

5.1. m = 3 (Ising model) 

The ground state averages {gk,i) show a different behaviour for k and / odd or even 
respectively. The values {gk,i) are positive for I — k odd (with a different behaviour for 
k odd or even) and zero otherwise. Thus for the continuum limit of the correlation 
functions one has to make a distinction between even and odd sites. From the analytic 
result for the semi-infinite chain (2.4) we find that for both correlation functions 

idKi) = (73^)^(^/0 • (5-5) 

The bulk limit is given by A; — )• oo with I — k fixed. Thus j — ^ 1 and we have x = ^ from 
(5.f ), corresponding to the Majorana fields. Note that the form (5.5) of the correlation 
functions is valid for arbitrary k and /. The scaling function F describing the infiuence 
of the boundary depends only on the ratio j. This can be expected of the continuum 
limit from general scaling arguments. In general the surface critical exponent Xs is given 
by [23] 

F{k/l) oc {k/lf'' for y ^ . (5.6) 

From (2.4) one has = | for the correlation function {g2j-i,2k) and = | for {g2j,2k-i)- 
Next we examine the numerical results extrapolated to the semi-infinite chain. The 
scaling behaviour of F can be seen from figure 1(a). In this figure we plot the contour 
lines {gk,i) {k — for k even and / odd. The straight contour lines passing through the 
origin imply a dependence of F on y only. The other correlation function with k odd and 
/ even shows the same behaviour. Assuming that F is a function of j we can calculate 
the bulk critical exponent x from the numerical data. Using extrapolated values for the 
functions {g2j,2k-i) and {g2j-i,2k) on the semi-infinite chain, we find x = 0.5000(1) in 
both cases. However, the surface exponent Xs can only be estimated with an accuracy 
of about 15% from the numerical data (extrapolated from chains of up to 40 sites), 
because the limit in (5.6) can not be reached with the available values k and /. 

5.2. m = 5 (three state Potts model) 

Next we consider the case q = exp(^), which corresponds to a three state Potts model 
as described in section 2. The computed {gk,i) are positive for I — k odd and negative 
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3 5 7 9 11 13 15 17 3 5 7 9 11 13 15 17 




2 3 4 5 6 7 8 9 10 11 12 13 14 



Figure 1. Contour lines of {gk,i)(l — k)^^ for the semi-infinite chain, (a) shows the 
data for q = exp(^) where x = ^, (&) the data for q = exp(i|^) with a; = |, and (c) 
the data for g = exp(i^) where x = — i. For (a) and (&) we use the values for k even 
and / odd with k < I < 18, for (c) we use all values with k < I < 14. The spacing 
between contour lines is 0.025 (a), 0.05 (h) and 0.01 (c). 
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otherwise. We find that we have to distinguish four non-zero correlation functions. That 
is, as in the case m = 3, the correlation functions described by \{gk,i)\ are different for 
k and / odd or even. The appearance of different functions for odd or even sites can 
be expected from the expressions for the correlation operators in terms of parafermions 
that arise from a quotient representation of the Temperley-Lieb algebra for q = exp(^) 
(section 2). 

From (5.2) we expect the bulk critical exponent x = |. In figure 1(6) we give the 
contour lines of F[k^ I) = g^^i (/ — h^^^ for k even and / odd. This plot is representative 
of the following: For the values A;, / < 18 (which we can determine the correlations of 
the semi-infinite chain for) none of the four scaling functions F can be described by a 
function of j. Furthermore we find that there is no value of x such that F is a function 
of J only. We conclude, that the range of A;, / is too small to be in the scaling regime. In 
summary, the critical exponents x and x^ cannot be computed in the case m = 5. This 
does not result from inaccuracies of computation and extrapolation of the correlations, 
but from the fact that in this case the chain length is not sufficient to give the continuum 
limit. 

5.3. m = I (Lee-Yang edge singularity) 

Here q = exp[^) and the hamiltonian H is related to a conformal field theory with 
central charge c = This non-unitary minimal model describes the Lee- Yang 

singularity of an Ising model in an imaginary longitudinal magnetic field [24]. Besides 
the identity it has only one primary operator with conformal dimension h = This 
is the dimension we expect as the bulk exponent from equation (5.2). 

The correlations {gi^m) show a behaviour different from the minimal models m = 3 
and m = 5. We find only one (positive) correlation function independent of A;, / even or 
odd. In figure 1(c) we plot contour lines for all values of A;, / < 14. The scaling function 
F[k^ I) = {gk,i){l — A;)~^/^ becomes constant for A; ^ 1 and / — A; ^ 1. Thus we expect 
F[k/l) = const, in the continuum limit yielding Xs = x. 

On the lattice there is an infiuence of the boundary, which vanishes quickly with 
k > 2. With the correlation function of the semi-infinite chain one obtains the bulk 
exponent from (5.1), taking A; > 3 and I — k > A and assuming a constant F. This way 
we find X = —0.20(1) verifying (5.4). 

6. Conclusion 

We have performed numerical investigations on the UqSU{2) scalar two-point operator 
for the UqSU{2) symmetric XX Z quantum chain at roots of unity. The properties of 
the correlation functions on the chain depend strongly on the value of q as we show for 
three different cases. For q = exp(^) one can find the (known) bulk exponent of the 
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two correlation functions {g2j,2k-i) and {g2j-i,2k) from small chains. For the non-unitary 
model given hjq = exp[^) it is possible to extract information about the continuum 
limit and critical exponents from chains with up to 24 sites. We have one correlation 
function and find the bulk critical exponent — | confirming x = hi^s = ^i^^- Thus 
the continuum limit of the correlations involves operators with spin hi ^3 as in the case 
q = exp(^), where this is already known from an analytical solution [2]. For the surface 
exponent we obtain = — | as well. For q = exp(^), where the XX Z chain corresponds 
to a three state Potts model, we find four different correlation functions {gk,i) depending 
on A;, / even or odd. In this case it is impossible to find the continuum limit from the 
lattice data (up to 28 sites). Thus the operators corresponding to the parafermions on 
the spin chain can not be found in this way. We have performed analogous computations 
for other integer m > 4. We do not present these results here, but want to remark that 
lattice effects dominate the correlations in these well. Thus in general it remains 

an open problem to find the conformal operators corresponding to the UqSU (2) invariant 
correlation operators for the quantum spin chain. 

Acknowledgments 

We would like to thank Vladimir Rittenberg and Fabian Ffiler for many helpful 
discussions. 

References 

[1] Pasquier V and Saleur H 1990 Nucl. Phys. B 330 523 

[2] Hinrichsen H, Martin P P, Rittenberg V and Scheunert M 1994 Nucl. Phys. B 415 533 

[3] Temperley H N V and Lieb E H 1971 Proc. R. Soc. A 322 251 

[4] Jimbo M, Miki K, Miwa T and Nakayashiki A 1992 Phys. Lett. A 168 256 

[5] Bonner J C and Fisher M E 1964 Phys. Rev. 135 A640 

[6] Borosowicz J, Kaplan T A and Horsch P 1985 Phys. Rev. B 31 1590 

[7] Gagliano E R, Dagotto E, Moreo A and Alcaraz F C 1986 Phys. Rev. B 34 1677 

[8] Grieger W 1984 Phys. Rev. B 30 344 

[9] Izergin A G and Korepin V E 1985 Comm. Math. Phys. 99 271 

[10] Karbach M, Mutter K H and Schmidt M 1994 Phys. Rev. B 50 9281 

[11] Batchelor M T and Kuniba A 1991 J. Phys. A: Math. Gen. 24 2599 

[12] Earner C J, Quispel G R W and Batchelor M T 1987 J. Phys. A: Math. Gen. 20 5677 

[13] Alcaraz F C, Barber M N, Batchelor M T, Baxter R J and Quispel G R W 1987 J. Phys. A: Math. 
Gen. 20 6397 

[14] Alcaraz F C, Baake M, Grimm U and Rittenberg V 1989 J. Phys. A: Math. Gen. 22 L5 

[15] Martin P P 1990 J. Phys. A: Math. Gen. 23 7 

[16] Martin P P 1991 Potts Models and Related Problems m Statistical Mechanics (World Scientific, 
Singapore) 

[17] Mittag L and Stephen M J 1971 J. Math. Phys. 12 441 



13 



[18] Yang C N and Yang C P 1966 Phys. Rev. 150 321 

[19] Bulirsch R and Stoer J 1964 Numer. Math. 6 413 

[20] Henkel M and Schiitz G 1988 J. Phys. A: Math. Gen. 21 2617 

[21] Belavin A A, Polyakov A M and Zamolodchikov A B 1983 Nucl. Phys. B241 333 

[22] Symanzik K 1981 Nucl. Phys. B190 1 

[23] Cardy J L 1984 Nucl. Phys. B240 514 

[24] Cardy J L 1985 Phys. Rev. Lett. 54 1354 

Appendix 

We list the ground state expectation values {gk,i) of the two-point scalar operators 
(2.1) extrapolated to the semi-infinite chain for the cases q = exp[^) corresponding to 
c = — Y (table 1) and q = exp(^) corresponding to c = | (table 2). Estimating the 
error as explained in section 4 these results are reliable up to four digits as given in the 
tables. We do not give the correlations in the case q = exp(ij^). They can be calculated 
from (2.4). 

Table 1. Correlation functions {gk,i) of the semi-infinite chain with A; < / < 14 for 
m = I (Lee- Yang edge singularity). The values given are extrapolated from chains 
with up to 24 sites. 

2 3 4 5 6 7 8 9 10 11 12 13 14 



1 1.81 2.392 2.788 3.108 3.384 3.629 

2 2.026 2.714 3.2 3.592 3.928 

3 2.041 2.737 3.23 3.627 

4 2.042 2.739 3.232 

5 2.043 2.739 

6 2.043 
7 

8 
9 



10 
11 
12 
13 



3.85 


4.054 


4.244 


4.418 


4 


578 


4 


746 


4 


879 


4.225 


4.494 


4.74 


4.965 


5 


177 


5 


381 


5 


551 


3.968 


4.269 


4.541 


4.788 


5 


015 


5 


236 


5 


42 


3.63 


3.971 


4.272 


4.543 


4 


793 


5 


024 


5 


224 


3.232 


3.63 


3.971 


4.272 


4 


544 


4 


793 


5 


012 


2.739 


3.232 


3.631 


3.971 


4 


272 


4 


544 


4 


784 


2.043 


2.739 


3.232 


3.63 


3 


971 


4 


272 


4 


544 




2.043 


2.739 


3.232 


3 


63 


3 


97 


4 


27 






2.043 


2.739 


3 


232 


3 


63 


3 


971 








2.043 


2 


739 


3 


232 


3 


63 










2 


043 


2 


739 


3 


231 














2 


043 


2 


734 


















2 


043 
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Table 2. Correlation functions {gk,i) of the semi-infinite chain with A; < / < 18 for 
m = 5 (three state Potts model). The values given are extrapolated from chains with 
up to 28 sites. 



k \ I 


2 3 


4 


5 


6 7 


8 


9 


10 


1 


0.9907 -0.1764 


0.3546 


-0.1064 


0.2139 -0.07631 


0.1514 


-0.05946 


0.1162 


2 


0.5136 


-0.1458 


0.1247 


-0.08158 0.06304 


-0.05589 


0.04036 


-0.04212 


Q 
O 




u .o^oy 


n 1 P.A7 




n 1 Q1 R 


— u .u ( oo 


n 1 '^Q7 
u . loy ( 


4 






0.584 


-0.1542 0.156 


-0.08962 


0.0822 


-0.06306 


5 








0.8051 -0.162 


0.287 


-0.09839 


0.1798 


6 








0.6139 


-0.1566 


0.1714 


-0.09219 


7 










0.783 


-0.1609 


0.2755 


8 












0.6309 


-0.1575 


9 














0.7694 



k\l 


11 


12 


13 


14 


15 


16 


17 


18 


1 


-0.04866 





.09376 


-0 


.04115 





.07825 


-0 


.03561 





.06696 


-0 


.03138 





.05834 


2 


0.02907 


-0 


.03358 





.02246 


-0 


.02778 





.01819 


-0 


.02361 





.01522 


-0 


.02047 


3 


-0.05799 





.1096 


-0 


.04806 





.09002 


-0 


.04106 





.07617 


-0 


.03586 





.06595 


4 


0.05365 


-0 


.04847 





.03902 


-0 


.03923 





.03028 


-0 


.03285 





.02456 


-0 


.02819 


5 


-0.0719 





.1319 


-0 


.05707 





.1044 


-0 


.04748 





.08643 


-0 


.04073 





.07369 


6 


0.0925 


-0 


.06561 





.06122 


-0 


.05092 





.0449 


-0 


.04155 





.03504 


-0 


.03504 


7 


-0.09741 





.1722 


-0 


.07113 





.1267 


-0 


.05649 





.1006 


-0 


.04706 





.08364 


8 


0.1808 


-0 


.09336 





.09908 


-0 


.06684 





.06623 


-0 


.05215 





.04891 


-0 


.04276 


9 


-0.1603 





.2679 


-0 


.09687 





.167 


-0 


.07072 





.1228 


-0 


.05611 





.0977 


10 


0.642 


-0 


.158 





.1872 


-0 


.09399 





.1037 


-0 


.06753 





.06984 


-0 


.05287 


11 







.7602 


-0 


.16 





.2626 


-0 


.09654 





.1632 


-0 


.07037 





.1199 


12 











.6499 


-0 


.1583 





.1918 


-0 


.09437 





.1071 


-0 


.06796 


13 















.7534 


-0 


.1598 





.2586 


-0 


.09632 





.1603 


14 



















.6557 


-0 


.1585 





.1954 


-0 


.09462 


15 























.7482 


-0 


.1597 





.2555 


16 



























.6604 


-0 


.1587 


17 































.7441 



